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Abstract 

We consider solutions of the 2-d compressible (isentropic) Euler equa- 
tions that are steady and self-similar. They arise naturally at interaction 
points in genuinely multi-dimensional flow. We characterize the possible 
solutions in the class of flows L°°-close to a constant supersonic back- 
ground. As a special case we prove that solutions of 1-d Riemann prob- 
lems are unique in the class of small L°° functions. We also show that 
solutions of the backward-in-time Riemann problem are necessarily BV. 



1 Introduction 

We consider systems of hyperbolic conservation laws in two dimensions: 

Most important are the 2-d compressible Euler equations for motion of inviscid 
fluids: U — {p, pv^ , pv^) {p density, v velocity) with fluxes 



f'={U)=v'=U + 



fy{u) = vyu + 



where p = p{p) is pressure. 

Our aim is to increase understanding of genuinely multi-dimensional flow, 
in particular its wave interactions. Some examples are regular reflection (four 
shock waves meeting at a point) l8l[17l[3|42l[231[15l[T3l[14]or Mach reflection 
(three shocks meeting with a contact or another type of wave) [Ij [2l [24l [25l 
[39l[37]. In these cases there are distinguished points near which the flow is, to 
leading order, constant along rays starting in the point. This leads to solutions 
that are self-similar and steady (from the point of view of an observer moving 
with the interaction point): 

U{t,x,v) = U{^), = Z(x,2/) e [0,2^[. 



*This material is based upon work supported by the National Science Foundation under 
Grant No. NSF DMS-0907074 
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In the case of Mach reflection the precise nature of the interaction remains 
controversial after decades of research. It is known that triple points (three 
shocks, with smooth flow in between) are not possible in most reasonable models 
(see [3T], [SI Section 129], [H], [3S1 Theorem 2.3]). However, beyond results for 
triple points and other special cases, the possible combinations have apparently 
never been classified systematically. Such a classification is our ultimate goal. 

We are particularly motivated by an example in |12| which features a steady 
and self-similar solution where two shocks and two contacts meet in a point. 
Numerical calculations suggest there is a second unsteady solution, with the 
steady one as initial data, so that the Cauchy problem for the 2-d Euler equa- 
tions would not have uniqueness, at least in its current formulation. Naturally 
we wonder which other steady self-similar solutions exhibit this behaviour and 
what characterizes them. 

The literature on multi- dimensional Riemann problems [291 I40[ 1411 128] is 
somewhat related to our flow class. However, in those problems only the initial 
data is necessarily constant along rays; we are interested in the special case 
where the forward-in-time solution equals the initial data. On the other hand, 
much of that literature focuses on initial data constant in each quadrant, a 
setting that is apparently so restrictive that the numerical studies have not 
encountered non- uniqueness phenomena like those observed in [T2]. 

In this article we focus on the case where C/ is a small (in L°°) perturbation 
of a constant supersonic background state U. Interestingly we do not need to 
assume that U is in BV, the space of functions of bounded variation; instead we 
will prove it (under standard assumptions about p{p))- This is crucial because 
in several space dimensions BV is probably too narrow to contain all reasonable 
flows [33], in contrast to one space dimension where a satisfactory theory has 
been based on BV or closely related classes (THl UHl 13] ■ 

Our results also apply to the classical case of 1-d Riemann problems for 
strictly hyperbolic conservation laws whose eigenvalues are either genuinely non- 
linear or linearly degenerate: for sufficiently small jump, their self-similar for- 
ward solutions (see [57] or [TSl Chapter 11] for construction) are unique in the 
class of L°° (rather than BV) solutions with small norm (for related uniqueness 
results see [11 Section 9.1] and [nHSlli 111301133 Hi I3H])- This generalizes 
an earlier result of Heibig |21| which required all eigenvalues to be genuinely 
nonlinear. While uniqueness need not hold backward in time, we are able to 
show that small-L°° solutions must be small-BV (which cannot be improved to 
any smaller commonly used class since examples with infinitely many jumps are 
easy to construct). 

Wc now summarize our main result. Consider the 2-d compressible iscntropic 
Euler equations. Let U £ L°° be a steady, self-similar, entropy-admissible weak 
solution, with IJJ/ — J7]]l°° < e for some supersonic background state U and 
e > 0. If e is suflaciently small, then U £ BV and it must have the structure 
shown in Figure [1] 
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Figure 1: U must be constant outside narrow sectors specified by eigenvalues 
evaluated at U. Linearly degenerate sectors: at most one contact discontinu- 
ity. Genuinely nonlinear forward sectors: at most one shock or simple wave. 
Genuinely nonlinear backward sectors: infinitely many waves possible, but no 
consecutive simple waves. Here we have taken the background state to have 
horizontal velocity (u, 0) and sound speed c. 



2 Balance laws 

Let V C IR™ be an open set. Consider smooth functions rj, ip^,ipy : V ^ IR. For 
^ C IR^ we say C/ = (;7\ ...,[/") G Ll^^{\R^ ; V) is a weak solution of 

v{U)t + i^^{U)^ + ipyiU)y <0 in A (1) 

if the inequality is satisfied in the weak sense (or: distributional sense): every 
X € A has an open neighbourhood N so that for nonnegative smooth $ with 
supp <I> g TV, 

- / <^tr]{U) + + %r{U)d{x, y,t)<0 (2) 

We call U a strong solution (or classical solution) if, in addition, it is a.e. equal 
to a Lipschitz-continuous function. 

Weak solutions — as well as other concepts — for the system of conservation 
laws 

Ut + r{U), + fy{U)y^O in A, (3) 

with f^, fy : V ^ IR™ smooth, are defined by interpreting ([3]) as 2m inequalities 
of the form (P), with = replaced by < or > and with ri{U) := U°', ip'^iU) := 
/"^"(C/), ^y{U) := P"iU) for a = 1,...,to. 
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We call [rj^ip^ ,Tpy) an entropy-flux pair for ([3]) if 

i'u=Vufu. ru=Vufl onV. (4) 

A weak solution J7 of ([3]) is called entropy solution (or admissible) if it satisfies 
([1} for all entropy-flux pairs with convex rj. However, all results in our paper 
hold even if we require ([IJ only for a single entropy-flux pair with uniformly 
convex 77. 

(An alternative to entropy-based admissibility are Lax conditions; see Sec- 
tion 15, especially (IHj) and (|45p . for Lax-type conditions suitable for our L°° 
setting.) 



3 Steady and self-similar solutions 

We arc interested in steady solutions: U is (a versio of a function that is) 
constant in t. Integrating by parts in t eliminates the first term in the integrand 
in <^ , after which using compact-in-t support and integrating with respect to t 
yields the equivalent statement 

^xr{U) + ^y^y{U)d{x,y)<0 (5) 



IR2 

for all nonnegative smooth compactly supported (now taken to be f-independent) 
functions $ : IR^ ^ IR. 

In addition we require U to be self- similar: a.e. equal to a function that is 
constant on each ray {s ■ {x,y) . s > 0} (for {x,y) ^ 0). While the derivation 
of the weak form is elementary, we show it in detail in order to point out some 
subtleties: To derive the weak form, first consider nonnegative smooth $ whose 
compact support is contained in the halfplanc {x > 0}. We may change variables 
in dSl) to (x,^) with ^ = y/x: 

0> - r [ [^A^,xOr{U{0)+%i^:^Or{U{0))xd^dx. (6) 

IR 

We take 

<l>{0 I ■ <^{x,x£,)dx = - x{<S>x{x,x£,) -\-^<S>y{x,x£,))dx 
Jo Jo 

so that 

hiO = / x<^y{x,xS,)dx, 
Jo 

oo 

x<i>x{x,xOdx = ~cl,{0~^MO- (7) 





^i.e. almost everywhere equal to 
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Then © is equivalent to 



IR 

for every smooth compactly supported nonnegative : IR — >■ IR, since every 
arises from ([7]) via 

where 77 is any smooth function with support in (0, 00) and integral 1. (|H]) is the 
weak formulation of 

{^yiu)-^r{u))^ + ^-^iu)<o. (9) 

By analogous calculations we obtain 

{Piu)'^nu))^ + r{u) = o. (10) 

If U is differentiable at ^, pO)) implies 

ifu{U)'m{U))U^ = 0- (11) 

If we repeat these arguments for x < 0, there is a single but crucial difference: 
the coordinate change to (|6]) produces an additional "— " from 

dy = \x\d£^ = ~x d^. 

The sign is irrelevant for ([TU]), but the entropy inequality ^ changes to 



{nu)-ir{u))^ + nu)\^o. (12) 



4 Smallness 

We restrict ourselves to the case where U is L°°-close to a constant background 
state U eV: _ 

\\U-U\\l^ < e. 

A finite number of times in this article, we choose e > sufficiently small for 
some purpose. 

4.1 Entropy gradient 

For any entropy-flux pair (77, -0^, ip^) and any w G IR™, 

fiiU) r]{U) + w ■ i^\U) := il}'\U) + w f {U) {i = x,y) 
defines another entropy-flux pair (r),'0^,'0^) since 

i^u ^^u + w ■ fh = ivu + w^)fli = Vufli {i = X, y). 



fjuu = Vuu, so convexity is not affected. By adding w'^ times (O to ([T]) (which 
is a hnear operation, hence compatible with weak formulation) we obtain 

fi{u)t + r{u)x + i'y(u)y<o 

which is equivalent to ([T} since we can reverse the operation using . Hence 
wc may assume, without loss of generality, that 

m(u) = (13) 

and do so from now on. 



5 Eigenvalues 

From (fTTj) . we see C/^ = is implied if the matrix 

is not singular. This suggests, as we show later, that U is constant in sectors 
where the matrix 

fuim-^iu) (14) 

is far from singular, so that the interesting behaviour is concentrated near ^ 
that satisfy 

= det {flj(U)-m{U)):=p{0 (15) 

for our constant background state U. The polynomial p has up to m real roots. 
Instead of focusing on one choice of coordinates, consider 

Pix:y)^det{xx fuiU)) (16) 

(with X = {x,y), / = (/^,/^) and (ai,6i) x (02,62) := aih - 0261) where 
we regard {x : y) <E IRIP^ as homogeneous coordinates; P is a homogeneous 
polynomial of degree < to. "x" is invariant under rotation, so that a coordinate 
change from {x,y) to [x' ,y') ~ R{x,y), R any rotation matrix, changes each 
root of P from {x : y) to {x' : y'). Each root ^ of p corresponds to a root (1 : C) 
of P. 

Since p has < m roots, wc can find some ^ which is not a root. The line 
(1 : ^) associated to ^ is, by rotating coordinates, aligned with (0 : 1). Then 

P(0 : 1) 7^ 0. (17) 

We assume from now on, without loss of generality, that this change has been 
made. 
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6 Change of variables 



6.1 Change to V 

(fT7|) also implies that is regular. Therefore is a diffeomorphism if we 
choose e > sufficiently small. Since we are using [/t = 0, it is not important 
to work with conserved quantities and we may change to 

V:=.r{U), V:=nU), 

and set 

fiV) P{U{V)). 

We let V be the open set of possible values for V from now on and abbreviate 
Ve:={veV \V~V\<ey 

pI7|) becomes 

if{V)~^V)i + V^O. (18) 
At points of differentiability ^ of ^ we have 

{fviVm-^l)VdO=0- (19) 

These are the same equations satisfied by a weak solution V oia. 1-d conservation 
law 

Vt + f{V), = 

that is self-similar, i.e. 

V{z,t) = ViO, 

if we identify ^ = z/t. Hence our x is a time-like variable while y is space- 
like. We could, for example, solve an initial-value problem by imposing data 
at a fixed x. However, there is no well-posedness without an entropy inequality 
which is what identifies the forward and backward directions of time in physics. 

6.2 Entropy-Flux pairs 

For any entropy-flux pair (77, t/;^, t/;^), define 

e{V):=riU{V)), q{V):^r{U{V)). 

Then 

ev = -ipuUv = VufuUv = VuVuUv = Vu- 

Therefore, 

evv = VuuUv = Vuuifu)'^- 



7 



We have 



qy = ipljUv = mffjUv ey/y, 



since 

fVl 



Therefore, {e,q) is an entropy-flux pah for (|f 8p . The entropy inequahty © for 
a; > becomes 

(<7(F) - ^e(F)) + e(F) < 0, (20) 



whereas 



for X < 0. 



{q{V)-^e{V))^+e{V)\^0 (21) 



6.3 Convex Entropy 

Lemma 1. If f^{U) is regular, then eyyr"r" ^ for all a = l, ..m. If f^{U) 
has only positive (negative) eigenvalues, then e is uniformly convex (concave). 

Proof. We shall use Proposition 6.1 from [36]. It states that if H is symmetric 
positive definite, and K is symmetric, then HK is diagonalizablc with real 
eigenvalues. Moreover, the number of positive (negative) eigenvalues of K equals 
the number of positive (negative) eigenvalues of HK. First, we write 

{iluu)~^ is symmetric positive definite, and eyv is symmetric. Then, apply- 
ing the proposition, since (/^)~^ is nondegenerate, all eigenvalues of eyv are 
nonzero. Moreover, if all the eigenvalues of {fij)~^ are positive (negative), then 
evv is positive (negative) definite, since a symmetric matrix is positive (nega- 
tive) definite if and only if its eigenvalues are all positive (negative). All that is 
left is to show that eyy?'"?'" 7^ in the case of eigenvalues of mixed signs. 

As in |341 Lemma 4.3.3], wc consider 

qv evfv- 

Then, 

qvv = evvfv + eyfvv- 

Therefore, 

evvfv ~ qvv — evfvv- 
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The first term on the right side is symmetric, and the second term on the right is 
a hnear combination of symmetric matrices, and is thus symmetric. Therefore, 
the left side is also symmetric and thus defines a symmetric bilinear form. Then 

evvifvrnr" = evvifw^y 

(Ac - Xp)evvr"rf^ = 0. 
Therefore, for (] ^ a, evyr"r^ = by strict hyperbolicity. Suppose that 

eyyr"r" = 0. 
By bilinearity, this would imply that 

eyyr^s = 

for all s € IR™. Therefore eyyr" must be the zero vector, but this contradicts 
the fact that eyv has all eigenvalues nonzero. Therefore, for each a, 

evvr°'r" 7^ 0. 

□ 



7 Versions 

Consider f - ^V{^) has a distributional derivative -~V € L°°, so 

there is a C £ IR™ so that 

fiViO) ~ ^V{0 t V{7j)dT^ for a.e. ^ & IR- (22) 

Jo 

Analogously, yields a C" G IR with 

q{V{0) - ^e{V{0) <C'- f e{V{v))dv for a.c. ^ e IR- (23) 

Jo 

Since the left-hand sides are continuous functions of V{S,) and the right-hand 
sides continuous functions of ^, Lemma [H] from the appendix applied to (|22p 
(with — split into <, >) and ((23)) yields a version (that is, an element of the L°° 
equivalence class containing V, which we will continue to refer to as V) of V 
that (a) has values in Ve everywhere, and (b) so that for all ^1,^2 G IR 

- V{r^)dr^ = (/(T/(6)) - 6^^(6)) - (/(^(eO) - ^iVi^i)) and 

(24) 

J^i 
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We abbreviate 

A{Vo + Ay, ^o) := / fv{Vo + s^V)ds (25) 
Jo 

and obtain 

/(Fo + AF) - /(K)) = A{Vo + Ay, yo)Ay, 

so 

^(6) - y(r7)dry = (a{V {^^) ,V - 6/) (y(6) - ^(Ci)) ■ (26) 

8 Strict hyperbolicity 

For the remainder of the paper we focus on the case of strict hyperbohcity. 
Many results would hold for weaker notions of hyperbolicity, but we prefer to 
keep the presentation simple. By strict hyperbolicity we mean that P in (jl6p 
has exactly m real roots {x : y) which arc necessarily distinct. That means 

det(/y(F)-C/) =0 

has m distinct real roots ^. 

A{V, V) = /y(y), so by smoothness of A we can take e > so small that for 
£ Pe there are m real eigenvalues X"{V^) {a ~ 1,...,to) of A{V^) which 
are smooth functions of and satisfy 

A"(y±) < A"+i(y±) vy=^,y± e T'e, a e {i,...,m- 1}. (27) 

(Vn is compact, so the separation is uniform, by continuity of A". The A" must 
remain distinct and real because their m real parts are continuous functions of 
V^, hence remain distinct for e > sufEciently small, so since A{V^) is real 
it cannot have non-real eigenvalues which come in conjugate pairs which would 
yield two equal real parts.) 

For a = 1, m wc choose a unit-length right eigenvector f-°'{V^) of A{V^) 
for eigenvalue A"(y^). f"(y^) is also a smooth function of V^. We choose left 
eigenvectors 1"{V^) that satisfy 

rf^=6a,p (a,/3 = l,...,m), 

which implies they are smooth as well. 
Abbreviate 

A{v) := Aiv,v) = fv{v), A"(y) :=A"(y,y), r"(y) :=r"(y,y), P(y) :=r(y,y). 
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9 Left and right sequences 

In this article we do not assume V G BV, so V need not have weU-defined left 
or right limits at any point Instead we consider pairs of sequences (C^T), (Cfe")i 
both converging to with ^^T < Cfc" ('^e do not require ^^T < ^ < ^j!" yet). Since 
V has values in the compact set Ve, there are subsequences {^^) of {^^) and 
(Cfc) of (Cfc) so that 

V{C)^V^ : Vi^k)^V'- (28) 
In such a context we write 

[g{V)] -.^ g{V+) ~ g{V-) 

for any function g (assuming there is no ambiguity as to which sequences are 
meant). 
Let 

Ji9{V);0 :=sup|[g(F)]| 

where the sup is over all sequences (^^) with the properties above. Then 
J{g{^):0 = if and only if 5 o y is continuous at 
By m. 



(i(F(e+), F(C)) - {V{^^) - ^(C)) = J^[^ Vi^+) - Virj)dr,. 
The limit as — >■ cx) is 

{AiV^)-Cl)[V]=0. (29) 
Hence for some a € {1, m} 

[V]\\f°'{V'^) and e = A"(F=^) (30) 
(that is, [V] is a scalar multiple of f^lV^ j). (P^j) is equivalent to 

[/(^)] - cm = (31) 

which is the usual Rankine-Hugoniot condition. Hence we may use any standard 
result that does not require continuity on each side of 5. 

10 General case 

In this section we collect results that do not require any assumption (such as 
strict hyperbolicity, admissibility, genuine nonlinearity, ...). 

Theorem 1. Suppose V is continuous on an interval I = ]^i,^2[ md that ^ is 
not an eigenvalue of A{V{^)) for any ^ G /. Then V is constant on I. 
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Proof. Fix some ^ € /. We claim that V must be Lipschitz at ^. Suppose not. 
Then we can choose a sequence {/i„} — > (with /i„ 7^ 0) such that 



< 



+ K) - V{0 



00. 



hn 

Divide both sides of ^ by + h„) -V{£,)\ to obtain 

+ K) - v{0 



1 



\V{£. + K)^V{i)\J^ 



\v{i + K)^vm 

V{i + K) - V{r])d-q = 0(1) as 71 ^ 00 (32) 



--0{hr,) 



{0{hn) since V is bounded). By assumption, A{y{£)) — is regular, so for h 
sufficiently small A(l/(^ + h), — S,I will be uniformly regular. That is. 



A{Vi^ + h),ViO)-il)v\>S\v\ 



3d>Q Vv e IR"^ 

Taking n — > 00, the left hand side of ([32]) stays bounded away from zero, while 
the right hand side goes to zero, leading to a contradiction. 

Therefore, V must be Lipschitz on /. Assuming ^ is a point of differentia- 
bility of V, we obtain 

However, as we assumed the matrix was regular on /, it follows that = a.e. 
on I. A Lipschitz function is the integral of its derivative, so V is constant on 
/. □ 

Theorem 2. Consider an interval I ~ ]Ci,^2[- There is a 6s — Ss{e) > 0, with 

6s iO as e 4- 0, 

so that 

Va e {1, m}Vx G / : |A"(V^(0) ^ £.\ > 6s (33) 

implies V is constant on I. [Here we do not require continuity of V , but a 
stronger bound on the spectrum.] 

Proof. Define 

6s -.^ sup |A"(l^)-A"(l^±)| 

and assume ([55)1 holds. The right-hand side converges to zero as e \ since 
A", A" are smooth and coincide for V = = V 

k 



Assume V is discontinuous at ^ G /. Then we may choose (^j!"), (^j, ) — J> ^ 
with l/(^^) -> and [V] ^ and obtain, by that 



e = A"(F±) 
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But then 

which contradicts 

Hence V is continuous on /; Theorem [T] yields the conclusion. □ 

11 Vertical axis neighbourhood 

As explained in the context of p7)) . we may choose some <^ G IR that is not a root 
oip (see (HH)) in the present coordinates and rotate coordinates so that (1 : ^) is 
aligned with (0 : 1) and therefore (0 : 1) with (-1 : ^). Then by is not 
a root oip in new coordinates, so Theorem [2] shows (if e > is sufficiently small) 
that U{t]) must be constant for t] in a neighbourhood of — ^. Rotating back to 
old coordinates it is constant — and therefore a weak solution — in sufficiently 
narrow open convex cones containing the positive and negative vertical axis. 
Therefore, we lost no generality by considering test functions supported away 
from the y-axis while deriving the weak form. 

12 Sectors 

By Theorem [51 wc can choose e so small that there arc intervals 

:=]A"(F)-(5",A"(F)+(5"[ (a = l,...,m) 

for 6" > so that V is constant outside U"=i-^"- We may choose 5" as 
e 4, 0. By forward sector (see Figure [1]) we mean ^ G /" with a; > 0, whereas 
backward sector refers to x < 0. 

13 Genuine nonUnearity 

Definition 2. We say /" is genuinely nonlinear if 

eV,:X'^{V)r°'{V)>0. (34) 

(if < we may without loss of generality flip the sign of r" (V) , f" (V^) (which 
remain unit-length) and 1°'{V)J°'{V^)). We say /" is linearly degenerate if 

eV,: A^(l^)r"(F) = 0. (35) 

14 Simple waves 

14.1 Simple wave curves 

Let s ^ i?"(y-,s) solve 

R"'{V-,0) = V-, i??(l/-,s) = r"(i?"(y-,s)). 

R" defines the a-simple wave curve. For each V~ we take the interval for s 
maximal so that R"{V~ , s) G Ve- 
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14.2 Wave fans 

If /" is genuinely nonlinear, then 

A"(i?"(F-,s)), = A?.(i?"(F-,s))r"(i?"(y-,s)) > 0, 

so 

s i-> \°'{R°'{V^ ,s)) is strictly increasing. (36) 
Let ^ 1-^ s(^) be its inverse map. By setting 

W{0:^R''{V',s{0) forC>A"(F-), 
we obtain a strong solution of ([T^ since 



=0 



If we interpret V as the value of W at the smallest ^, then only the s > part 
of R", denoted i?""*", is relevant. 



15 Discontinuities 

We recall some standard results we need later, to show that they do not depend 
on having a smooth neighbourhood on each side of a discontinuity. 

15.1 Shock curves 

Consider sequences (^^1") and (^^7) converging to ^, with < for all k, 
so that V{^f) -> V^. This is the setting of ^ which implies [V^] is a 
right eigenvector of A{V^) and ^ the corresponding eigenvalue. So there is an 
a € {1, m} with 

h{V+,s) := - - sf^iV^) = 0. 
/i is smooth, h{V~,0) ~ and 

dh 



so the implicit function theorem yields, after taking e > sufficiently small, 
existence of a smooth bijcctivc map s i— > S"(V~ , s) with 

s°'{v-,o) = V-, s°'{v-,s) -y- - sf'^iy^) = o. 

For each V~ we take the interval for s maximal so that S°'{V~,s) £ Ve- 
S°'{V~,-) defines the a-shock curve of . It contains (via s = 0) and 
has tangent r°'{V~) there. 

We take e > so small that for each a only = S°'{V~ , s) are solutions 
of dSIl). 
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15.2 Contact curves 

Assume /" is linearly degenerate. Then 

A"(i?"(F-,s)), =A^(i?"(l/-,s))r"(i?"(y-,s)) #0. 

Hence 

s A"(i?"(F-,s)) is constant. (37) 

Now consider 

F{s) := f{R"{V-,s)) - f{V-) - aR''iV-,s) - V"). 
Then F(0) = 0, and 

This is zero if we set ^ = A"(i?"(F^, s)) which is possible since the latter is 
constant. Hence the Rankinc-Hugoniot condition (PT|) is satisfied. 

Since i?" is maximal in V^, since 5" is maximal as well and contains the 
only points in satisfying pip , and since both arc simple smooth curves, they 
are identical. 

Hence, at ^ where an a-contact — [V] a right eigenvector for X°'{V^) — 
occurs, we have 

X°'{V-)^S, = X°'{V^)^\°'{V+). (38) 

15.3 Admissible shock curve 

Now assume /" is genuinely nonlinear. Assume V is admissible. Consider a 
forward sector first. The entropy inequality 

[q{V)] - ^[e{V)] < (39) 

can be derived from ([20]) in the same way as ([31]) from ([XSl - 

By ([30]), a jump from V~ to V+ must be located at ^ = A"(F+, V~), and 

X^(V~,V~) = A"(F~), 

so 

diX'^iV-, V-) + d2X"iV-,V-) = X'^iV-) (40) 

Moreover, 

X''{V-,V+) = x°'{v+,v-) 

since 

A{V+,V-)^ f fv{{l-s)V-+sV+)ds= f fv{rV- + {l-r)V+)dr ^ A{V-,V+). 
Jo Jo 
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Therefore 

di\'^{v-,v-) = d2\''{y-,v-). (41) 

Combining and (|1T|) we have 

9iA"(y-,y-) = ^^^^ -a2A"(v^-,v"-) 

Thus 

A"(5"(y-,s),F-)3 =°aiA"(v^-,v^-)5ni/",o) = iA^(y-)r"(F-) >0. 

Hence for e > sufficiently small 

s K> A" (5" (F" , s) , F" ) is strictly increasing. (42) 

We may reparametrize the a-shock curve of to be A" = ^ i—> W{£). 

Abbreviate := \°'{V~). To avoid clutter we change coordinates so that 
e(y-) = 0, q{V-) = 0, f{V-) = 0, VF(Co) = ¥' =0 (which is acceptable since 
adding constants to V, /, e or q has no effect in and ([201) )• (|3ip becomes 

o = /(mO)-^w^(0, 

with derivative 

0- (/v(W^)-e/)We- PF. (43) 
((39)) is equivalent to E{C) < for 

We analyze the situation near ^ = ^o- Since 

Ei^o) = q{W{^o)) - ^oeiWi^o)) = q{0) - $06(0) = 0, 
so we need to consider the first derivative, given by 

= {qv - ^ev)W^ -e = ev{fv~ ^I)W^ - e # eyW - e. 

Then 

EdCo) = ev(M^(eo))M^(eo) - e{W{^o)) = 0. 

=0 =0 

Hence we need to consider the second derivative as well: 

^44 = evvW^W + eyW^ — eyW^ = eyyW^W. 

Then 

=0 
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The third derivative finally yields a result: 

E^^^ {evvvW^W^ + evvW^^)W + eyvW^W^ 

so 

Ee^ao) evv{W{^o))W^{^o)Wao) + 0, 

because W(^o) — V~ and W^^(Co) — r°'{V~) by definition of the shock curve. 
Hence, if eyyr^r" > 0, for e > sufficiently small, 

E <0 ^ ^ < Co- 

Therefore, in this case, only the s < part of the shock curve (corresponding 
to C < 1^0 due to A" strictly increasing) yields admissible shocks. We call this 
part 

If eyyr"r" < 0, then the £, > £,o part is relevant. 
To this end, if 

eyvr"r" > 0, 

we define the "forward sector" to have x > 0, and the "backward sector" to have 
a; < 0. Since everything is smooth and this quantity can never be zero, it must 
be positive for all V~ G Ve if it is positive anywhere (vice versa for negative). 
Conversely, if 

eyvf"r" < 0, 

the "forward sector" has a; < and the "backward sector" has x > 0. 

We can consider the same setting but for an a; < sector: an analogous 
argument, starting with the opposite entropy inequality (21). 

Moreover, (42) shows that admissible shocks in forward sectors satisfy the 
Lax condition 

x{v-) > e > A(y+). 

More precisely the following uniform Lax condition holds for forward sector 
shocks: there is a constant (5l > so that 

X{V-)~Sl\[V]\>^>X{V+) + Sl\[V]\. (44) 

Finally we consider the same setting but for a backward sector: 

AC^"") + '^l|['^^]| < ( < X{V+) - Sl\[V]\. (45) 

If the background state for Euler flow is supersonic horizontal velocity to 
the right; then all forward sectors are a; > 0. If, however, the background state 
is supersonic horizontal velocity to the left, then all forward sectors are a; < 0. 
For 1-dimensional conservation laws with convex entropy, all forward sectors are 
t > 0. Finally, it is not hard to construct examples that satisfy all assumptions 
but have possessing eigenvalues with different signs, so in certain cases there 
can be forward sectors for some eigenvalues contained in a: > and others in 
a; < 0. 
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16 Linearly degenerate sectors 



Wc consider linearly degenerate /" and allow both x > and a; < 0. 

Lemma 3. (a) A" o is continuous, (b) If S, ^ -^"(^(C)) on an open set 
A G I" , then V is constant on A. 

Proof. Assume o V and therefore V are discontinuous at £ I". Then we 
can choose (^^) with V{^^) so that 

[X'^iV)] ^ 0. 

However, since /" is linearly degenerate and is on the a-simple wave curve 
of dnS]) shows 

A"(F+) = A"(F-)=eo, 

contradicting the assumption that A" o F is discontinuous at .^o- This shows (a). 
Theorem [1] yields (b). □ 

Lemma 4. For any subset E d I"' , and for almost every € there exists 
D <Z E containing such that 

(^r>)'(Co) exists and is finite. 

Proof. The idea is to use [161 Corollary 1] to obtain differentiability after restric- 
tion to a subsequence. However, the result on which the Corollary depends is 
only true for functions from a subset of IR" to IR™, with n > m, which need not 
be the case for our V : \R ^ IR™. We instead apply the result to the following 
function from C IR to IR. For the background state V and ^ G define the 
function 

[IB] proves that for any E C I", for almost all G E there exists D' with 
EZ) D' 3 ^0 such that 

(P(y)V|£)/)'(^o) exists and is finite. 
Recalling ([26]). we have 

For ^ 7^ a multiply 1^ (y{^),V{^o)) on the left to obtain 

•'Co 



We then estimate, for M only depending on V^: 

A/|c-eol > \\^{V,V^)~^f{V,Vo){V{0~V{io))\ 
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(HZD bounds |A'^(y, Vo) - Col from 0; so, with some other constant M' we have 
TliGrGforc 

^^i^{v{0,v{^o)){v{0-v{^o)) 

(which is at C = Co) is Lipschitz at Co with constant < A/'. This implies that 
the difference quotients 

i^{v{0,vi^o)){viO - Vj^o)) - ^^(^(eo),^(Co))(^(Co) - V(Co)) 

C-Co 

are contained in Bm'{0), a compact set in IR™ (for C sufficiently close to Co)- 
Therefore, for each 13 ^ a we can successively pass to nested subsequences in 
D' so that we finally obtain D such that 

{r(V)V\D)\^o), {i^{V, Vi^o))V\D)'iCo) exist and are finite 

(where D C D' C E (Z I" and S,o e D). We now claim that 

W^g^iW) ■.= 1"'{V)W 

and 

W ^ g^iW) := i^{W, F(Co)) {W - ^(Co)) (/? ^ a) 

yield a local diffeomorphism V, 3 W g{W) {g\W), g'^'iW)) (for 
possibly smaller e). To see this, notice that for /3 ^ a 

= 9^^{V{^o))z^l^{V{^o))z ^ z||r"(y(Co)). 

Then, 

= g^{V{^o))z^r{V)z = r{V{(o))z + 0{e)z z = 0. 

Since (g o V)\d is differentiable at Co and g is a local diffeomorphism, we have 
that is differentiable at Coi ^^nd the lemma is proved. □ 

Theorem 3. On a linearly degenerate (forward or backward) sector, V is either 
constant, or constant on each side of a single contact discontinuity. 

Proof By LemmaO F -.^ e I" | C = A"([/(C))} is closed and V is constant 
on 

Assume there are C11C2 G F and rj G I" with Ci < < C2- Then we can 
choose a maximal Jry^, rj^l containing r] but not meeting F. Necessarily 77^ £ F, 
so ?7+ = A"(y(?7+)) and 77^ = A"(y(77")). But V is constant on ]r]~,ri'^[, so 
77+ = rj^ , which is a contradiction. 
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Hence F must be a closed interval. 

Assume F has positive length. By p4| . f{V) — Sy is Lipschitz and therefore 
differentiable on. E C F (where F \ E has measure zero) with 

{fiV) - (V)^ + V^0 onE. (46) 

Note that this is in the strong sense, not just distributionally. By Lemma |31 
for almost every ^ G E wc can find D C E containing ^ such that V\d is 
differentiable at .J. Thus if F has positive length, there exists ^ & D C E C F 
such that (Uni) holds and V^jj is differentiable at ^. Therefore, we have 



sod^V\D{0 II r"(y(0), hence 
by linear degeneracy. However, 

e = A"(y(c)) 

implies 

l-A^(y(0)9?Mi5(?), 

which is a contradiction. 

Hence F must be a point (or empty, which can but need not be ruled out). 

□ 



17 Genuinely nonlinear sectors 

Consider a genuinely nonlinear /" . Consider either the forward or the backward 
sector. We partition /" into the three sets 



s 


= {e e 


J(F;e) >0}, 


n 




j{v;o^o. e-Afc(F(e))}, 


c 







where S stands for "shock" , TZ for "resonance" , C for "constant" . Complements 
(denoted by C) are taken with respect to 



17.1 Backward sectors 

Consider a backward sector (sgn(a;) = — sgn(ev'v'r"r")). Assume V is admissi- 
ble. 

First we observe crucially that shocks of admissible V must have a left and 
right neighbourhood in each of which V is constant. The neighbourhood size is 
lower-bounded proportionally to the shock strength. 
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Figure 2: For admissible V ^ each shock has a constant neighbourhood with 
lower size bound proportional to the shock strength. Reason: only in-admissible 
shocks could jump A(y(^)) back to ^ immediately. For backward sectors, con- 
secutive shocks or shocks interspersed with compression waves are possible. 

Theorem 4. For any Cz S there are cr+(^o) > Co (maximal) and cr^(^o) < S,o 

(minimal) so that V is constant on [c^ {S,o),£,o[,]£,o,cr'^{£,o)] C Moreover 
a±(Co)e7^ua/", 

<J-{^o)<Co-SLJ{V;io), (50) 

cT+{Ca)>(o + SLJiV;^o), (51) 

and 

^"{V{(o+))-^o>SlJ{V;^o), (52) 

>^"{Vi^o-)) - eo < -<5l J(F; ^o) (53) 

(where Sl is as in (j45p ). 

Proo/. (See Figure H) 

1. Assume V is discontinuous at ^o- Then we can choose a strictly decreasing 
sequence {(,'^) I and another sequence (^^) — > so that < and 
V{^^) The backward Lax condition ^ implies X°'{V+) - > 0. 

Assume there is no 5 > so that - ^ > for ^ G ]^o, Co + S[. Then 

we can rename (^^) to (^^T) and V'^ to (replacing the previous choice) 
and choose a new decreasing sequence {^^) I so that A"(y(C^!")) — Ca!^ < 
and so that V{^'^) V~^. We may assume, by omitting members from both 
sequences, that < for all k. Then the backward Lax condition ps)) yields 
A"(y+) - Co > 0, but that implies A"(y(C^!")) - > for /c sufHciently large, 
which is a contradiction. 

Thus we may choose a maximal (7+ (Co) <= -f" H ]Co, oo[ so that 

VCe]eo,^+(Co)[: A"(l^(e))-C>0. (54) 
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Analogously we obtain a minimal a (^o) <= -f" H ] — oo, so that 

vce]a'(Co),eo[: A"(i^(e))-e<o. (55) 

2. If ]^o, ""^(^o)] contained a ^ £ 5, then we could choose 

'?€]eo,^+(eo)[n]a-(0,e[ 

so that 

mil for ia l55ll for ? 

A«(F(ry)) > > A"(y(ry)), 

which is a contradiction. Hence 

y is continuous at every ^ e ]^o, "'''^(Co)]- (56) 

By Theorem [U ([5^ combined with ^ 7^ '^^(^(0) for /3 7^ a (by definition of 
J") yields 

V is constant on ]^o, cr'''(^o)]- (57) 

Analogously we show V is constant on [o"^('^o), Co[- Then we may take any 
(^fe ) Co with <i^< i'l and and V{if)^ and obtain (gl]) and ([53]) 
from the backward Lax condition p5|) . 

3. The boundary cr'*'(Co) with the property (|54l) is maximal. If it is not a 
boundary point of then there is a sequence (r/„) 4, (t+(^o) in with 

A"(^(?/n))-?7n<0. 

By ([55)) that means 

A"(V^(a+(eo)))-cT+(eo) <0. 
On the other hand, ([5^ and ([55]) show < is not possible, so 

a+(^o) = A"(y(a+(eo))) 
which means (t"''(^o) G ^-nd 

1^ 113 

a+(eo) = A"(y(a+(^o))) = A"(nCo+)) > eo + '5LJ(i^;e) 

which implies (|5T]) . Analogously we obtain ((50)) . □ 

Remark 5. In particular 5 is discrete, hence countable. (This does not imply 
V G BV(/") yet until we also show the continuous part of V has finite variation.) 

Since we have shown now that V has well-defined left and right limits in 
each discontinuity, we may modify V in each ^ G 5 to be the right limit. 
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Lemma 6. There is a constant Cg, independent ofV, so that for any G S, 

C^]a-(Co),a+(Co)[ 

implies 

J{V;^o), |A"(l/(^o+)) " Co|, |A"(V(eo-)) - 6| < Cs\i - U (58) 
Proof. (Sec Figured) 

tsot 

l^-^ol >min{|a-(^o)-?o|,k+(eo)-eo|} > '^lJ(F;^o) 

^ ^(V;eo)<'5Z'ie-eo| 

|A"(l/(eo+))-$oUA"(y(eo-))-eo| <'5LJ(y;^o)< ie-eo|. 
Take C5 > max(l,(5^^). □ 
Lemma 7. If ^ is a limit point of S , then ^ G 72,. 

Proof. (See Figure O) Let (^„) — ^ be a strictly decreasing sequence in 5 (the 
strictly increasing case is analogous). S is discrete, so ^ ^ [^n) i {£,n)[ (it 
could not be a limit point otherwise). Choose some ?7„ £ ]o'~(Cn)jCn[ each 
n. Then 

|A"(F(r?„)) - Vn\ = |A"(F(C„-)) - r^l (59) 

< |A"(nC«-))-en| + len-^n| 

< c5ie„-ei + i^n-ei™o. 

(''/n) C and ^ ^ 5, so A"ol/ is continuous at ^ and therefore A"(F(^)) = ^. □ 

Theorem 5. //^o G C, t/ien V is constant on an interval ]k~ (^q), k,'^{S,o)[ that 
contains ^o- We take the interval maximal in k^(^o) either in TZ Li S or 
endpoints of I" . 

Proof. By Lemma [71 ^0 is not a limit point of S (since it would be in TZ other- 
wise, and TZO C = 0. Hence V is continuous in a neighbourhood of ^o- Then 
A"(V(Co)) - Co ^ implies A"(F(0) - e ^ for C in a neighbourhood of fo- 
Since A'^(F(C)) - ^ 7^ for ^ ^ a by definition of Theorem [T] shows V is 
constant on this neighbourhood. We may take ]k~(Co), «:^(Co)[ as described in 
the statement. 

By what we have already shown, k^{£,o) ^ C because it would violate their 
extremality. □ 

From now until just before Theorem [71 we consider only V restricted to 
without writing V\ja to avoid clutter. 

Lemma 8. // G T^, then X" o V^q^ is Lipschitz at with Lipschitz constant 
<Cs + 2. 
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Figure 3: ^ i->- X{V{^)) — £, is Lipschitz at points where it is 0, since shocks 
have to weaken at least proportionany to their distance from ^o- 



Proof. (See Figure[3l) Consider £, ^ S with ^ > (the case ^ < is analogous). 
We first obtain a Lipschitz estimate for A"(V^(^)) — ^. 

1. If ^ e 7?., then by definition of TZ 

|A"(ne))-eho. (60) 

2. If ^ € C, then k~{^) G TZUS by Tlieorem[5l cannot be a boundary 
point of /" since < < CO 

2a. If n-{C) e n, then \°' {V (k^ {^))) = n-{^), so 

A"(F(0) = X" iO)) - ^-{0 e [eo,e] 

^|A"(ne)) - - el < iCo - ci- (ei) 

2b. If K-{0 e 5, then (note ^ ^ C since e 

7?. n C = 0, so Lemma [5] applies to ^o) 

|A"(F(c))-e| = |A"(F(«-(e)+))-e| 
< |A"(F(^-(c)+)) - K-(e)| + i«;-(o - CI 

Lemma [6] 

< C5|a^-(0 - Col + k-(0- CI <(^^5 + 1)1^0 -C|. (62) 
Combining all cases (|61l) . (15^ we see that 

|A"(y(c)) - e - (A"(^(Co)) - Co) h lA"(^^(e)) - ci < (c^ + i)ieo - ci. 

V ' 

=0 

Hence A" o V^^^ itself is also Lipschitz at with constant < Cs + 2. □ 
Theorem 6. V is Lipschitz with constant < C-jz (independent ofV) at any 
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Proof. For each p ^ a multiply £^ onto to obtain, for some M depending 
only on P^. 



M\i-U > ^^^{V{0,V{(,o))[A{V{0,V{i^))-S,^l){V{0-V{^o)) 

= (A'5(y(0,^(eo))-A"(y(Co)))^"''(l^(0,n6))(^(C)-neo)) 

Since A^ — A" is bounded away from by ((27)) . we obtain for some other constant 
M' that 



so that 



(which is = at ^ = 6) is Lipschitz at ^ with constant < M . 



W ^ g^{W) := £^{W, V{Co)) {W - ^(^0)) (/? ^ «) 



and 



g°'{W) A"(M^) 

yield a local diffeomorphism V, 3 W ^ g{W) := {g^(W), g'^iW)) (after 
reducing e > 0, if necessary): for /3 7^ a, 

= 3^(n^o))^ = ^^(n^o))^ 
implies z || r"{W), and then 

0-.9S^(^(eo))^ = A^(F(eo))^ 1* z = (63) 

by genuine nonlinearity; hence 31^(1^(^0)) is regular. 

Lipschitz continuity at ^ for C5 9 ^ i-> g{V{^)) (from Lemma [8] for g") 
implies Lipschitz continuity at ^ for C5 9 ^ i—> V{^). □ 



Lemma 9. Define 



(64) 



Then Vs is well-defined and a right- continuous saltus function (see Definition 
\13\ in the Appendix). 



Proof. 



(Uni , ^ 

< (2<5i)-i E - ^ (2'5l)-'|/"I < 00 
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since the neighbourhoods ] cr~ (77), tT+ (77) [ oi i] E S are pairwise disjoint and 
contained in Hence not only is S countable, but the jumps sum to a finite 
number. Hence (|64p makes sense. In Definition [T3] only &„ are used so that Vs 
is right-continuous. □ 

Lemma 10. Vs is Lipschitz with constant < Cs (Cs independent ofV) at any 

Co 

Proof. Let ^ > ^0 (the case ^ < ^0 is analogous). 

1. If ^ ^ ]a~ {rj) , {r])[ for some 77 € 5, then we may estimate 

< {26l)-'\^ - (65) 

since the neighbourhoods ]a~ (rj), {r])[ of distinct r] € S are pairwise disjoint 
and contained in [^,^o]- 

2. If ^ G ]<T^(r/), cr+(r/)[ for some 1] E S, then we apply (|65|) with ^ ^ a^{rj): 

\Vs{0 - Vs{£.o)\ < \VsiO - Vs{<j~{rm + \Vs{'y-{v)) - Vs{^,)\ 

< 1^5(0 - ^5(^^(^))| + {26L)-'\<y-{n) - Col 

< 1^5(0 - Vs{a-m + {2Sl)-'\^ - U 

2a. For ^ e ](j^{ri), r/[ the first term is = 0. 

2b. For ^ e [7],a^{7])[ the first term is estimated by Lemma |6] (using ^0 ^ 
]a-{r]),a+{r])[ C C since ^0 e 7^, 7^nC = 0): 

|^s(0-^^5(^-(r?))| - J(^^;^) < C^5|^-eo|<C5|^-Co|. 

Altogether we get the desired estimate, with Cs := C5 + (25^)^^. □ 

Theorem 7. V = Vs + where Vl is Lipschitz, with a Lipschitz constant 
independent of V . In particular V is BV. 

Proof. It is sufficient to obtain a Lipschitz estimate for ^, 77 G /" (^ < 77) since V 
is constant in between intervals I^ , the distance to for (i ^ a has a positive 
lower bound independent of V , and V is bounded. 

1. First consider ^ G TZ. 

\VLiv) - VLiOl < l^(^) - viO\ + \Vsiv) - Vsm 

Theorem [gI 

< C|77-^| (66) 

Lemma llOl 

for some constant C independent of V. 

2. Now consider ^ G S. Then is constant on ]o'^(^),(T+(^)[ (the jump of 
y at ^ is cancelled by Vs), so we only need a Lipschitz estimate for rj > o'+(^) 
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(which impUes cr~^{^) ^ 31°'). By Theorem HI cr^(C) S so we may use 
(with e ^ a+(C)) and Vl{£.) = YUa+iO) to get 

I^L(e) - = \VL{a+{0) - VlMI < C\a+{0 - rj\ < C\i - (67) 

3. FinaUy consider ^ e C. Then Vl (hke V) is constant on ]k^(^), k+(^)[, so 
we only need a Lipschitz estimate for rj > k+(^) (which imphcs k^{£,) ^ 81°'). 
By Theorem [5J k+{^) e 7^U 5. 

3a. For k+{^) € 7^ we may use §^ with ^ ^ k+{^) e TZ and ¥£,(0 = 
^^(^^+(0) to get 

|^l(0 ~ = |14(A^+(e)) - ^l(??)| < C|k+(C) - ryl < - (68) 

3b. For K+(0 e 5 we may use with £_ ^ k+{£_) e S and V^l(0 = 
yL(«;+(0) to get 

\VLiO - VlMI = \VLin+{0) - yL{^)\ < c\k+{0 - v\ < - vl m 

□ 

Remark 11. This shows that entropy- admissible self-similar weak solutions to 
the Riemann problem (for sufficiently small jump) for 1 — d hyperbolic con- 
servation laws are unique in L°° (assuming \\U{-) — C/||l°c sufficiently small), 
extending the well-known result that they are unique in BV (see Theorem 9.4.1 
in [in]). 

17.2 Continuity on open nonempty intervals 

Theorem 8. Consider any genuinely nonlinear sector, forward or backward. 
If V is continuous on an open interval B G 1° , then it is either constant or 
constant on either side of a single a-simple wave. 

Proof. (See Figure [SJ) V is Lipschitz on B, since we can repeat Lemma |S] and 
Theorem [6] with obvious changes to their proofs {S need not be considered since 
V is continuous here). 

By continuity of C H i? is open, hence a countable union of disjoint open 
intervals. V is constant on each of these intervals, by Theorem[Tl and so is X^oV, 
so that X°'{V{(,)) — ^ = cannot be satisfied at both endpoints. Therefore at 
least one cndpoint of each of these intervals is not in TZ. Suppose this endpoint 
is in iS. Since SClB is empty, this endpoint is an endpoint of i?. If this endpoint 
is not in 5, then it still must be an endpoint of B. Since there are only two 
endpoints, C H i? is a union of at most two of these intervals. 

It follows that C n B is either B itself, B minus a single point (which by 
continuity of V implies V is constant on all oi B), or B\{TZn B), where TZO B 
is a closed interval of positive length. By definition of TZ, 

A"(nO) = e (70) 
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onTZriB. By ([TO)) (defined a. e., since F is Lipschitz) is a multiple of 
Therefore ^ i— >■ V{£^) is part of the a-simple wave curve i?", and (|70p shows it is 
the ^-parametrization of R". Hence V is an a-simplc wave on TZ. □ 

Remark 12. This shows that although infinitely many waves can occur in a 
backward sector, there cannot be consecutive simple waves. For more than one 
simple wave to exist, there must be at least one shock in between. 

17.3 Admissible genuinely nonlinear forward sectors 









shock 







Figure 4: In a forward sector ^ i— >■ 
A(F(^)) — ^ cannot return to after 
a shock, and has the wrong sign for 
another admissible shock. 
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Figure 5: In a forward sector, after a 
simple wave ^ i— )- A(F(^))— ^ has the 
wrong sign for an admissible shock, 
so it cannot return to 0. 



Theorem 9. Consider an admissible genuinely nonlinear forward sector. Then 
V is either constant, or constant on either side of a single simple wave, or 
constant on either side of a single shock. 

Proof. (See Figures 0] and [51) Assume that V is discontinuous at some G 
Choose (^-), iCt) Co with < C and V{^^) where [V] ^ 0. The 

forward Lax condition pi)) yields 

A"(F-)>eo>A"(F+) (71) 

We may proceed in the same manner as in the proof of Theorem [H X{V{£,)) 
is still constant in ]Co,cr^(Co)[ and ^ is strictly increasing, but now ([71]) has 
the opposite comparisons: ^"(V {£,+)) — £, is negative and cannot reach or 
change signs again. Hence a-~^{^o) is the right boundary of By an analogous 
argument on the £, < £,o side we obtain that cr~(Co) is the left boundary of 

Now assume V is continuous on Then Theorem [8] yields the rest of the 
result. □ 
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18 Isentropic Euler 
18.1 Calculations 

We now focus on a particularly important case, the isentropic Euler equations: 

Ut + nu)., + fy{u)y = o, 



u = 



m \ n 



mnp ^ I \ n^p ^ + p 



Here {m,n) is the momentum density vector, v = = (^, ■^) the velocity. 

Then P is an open subset of {(p, m,n) g IR'^ : p > 0}. We assume the pressure 
P = p{p) satisfies 

c^^p'{p)>0 

for all (0 > 0; c is the sound speed. We assume 



c 



> -1 (72) 



which is satisfied for most relevant pressure laws, including = p'' for 7 > — 1. 
Take 



then 



P 



viU) := p{e{p) + l\v\'), m) ^ (rKC/) + p)v 



form an entropy-flux pair {ri,ip) with uniformly convex 77. 

For simplicity we assume units have been chosen so that c ~ I for p = 1. 

The Euler equations are invariant under rotation (and mirror reflection): if 
i7 is a weak/weak entropy /strong solution, then for any 2x2 orthogonal matrix 

U'^{p\v'), p'{t,x') ^ p{t,x), v'{t,x') = Qv{t,x), x' = Qx 

is another weak/ weak entropy/strong solution. The equation also also invariant 
under change of incrtial frame: for any a E IR^. another solution U" is 

U" = ip",v"), p"it,x") ^ p{t,x), v"{t,x") = v{t,x) + d, x"^x + ta. 

Consider steady self-similar solutions. In the framework of the present pa- 
per we consider only the strictly hyperbolic case. To this end we consider a 
background state U = {pq,Mo,0) with Mq > 1. (Due to rotation invariance no 
generality is lost. If we interpret supersonic steady Euler flow as an initial- value 
problem, with data imposed at x = —00, hyperbolicity with x as time and y as 
space variable requires M > 1, not just |M| > 1.) 
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In addition we choose e > so small that j|J7 — < e implies M > 1 as 
well. We may also choose units so that po = cq = I. 
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The eigenvalues of /^(l,Mo,0) are A/q ± l,A^o- Therefore, if A/q > 1 as as- 
sumed, all eigenvalues of will be positive, making ew positive definite, and 
therefore the forward sectors have a; > 0. If instead Mq < — 1 (as required for 
hyperbolicity with —x serving as a time variable), then all the forward sectors 
would have a; < 0. 

The generalized eigenvalues (roots of p in (fT5|) ) are 



Ah 



{pcf 



{pcy 



Ao 



n 

to' 



which are real, distinct and analytic functions of U for Af > 1. 
The generalized eigenvector vq for Aq is (0,TO,n). 



Vc7Ao(?7) • ro{U) = V(p,M,jv) (-) • (0, m, n) = (O, 



-n 1 



) • (0,TO,n) = 0, 



so the 0-field is linearly degenerate. 

For the zb-fields it is sufficient to consider the generalized eigenvectors only 
at U: 





" 10" 







1 " 




±TO 




1 - 2m 




—mn n 


TO 


- 

H- 


±(to2 - 1) 




—mn n m 




1-n^ 


2n 




V m"^ — 1 



dX± m± pcn{m^ +n^ - {pcf)-'^/'^ p=c=i,n=o M 



dn m^ — (pc)2 Af2 — 1 

For the p, m derivatives we may substitute n = first: 

A± = , = (73) 

v/m2 - (pc)2 VA/2 - 1 
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Then 
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Therefore, 



Vc/A±([/) • r±([/) = - 1)--V2 



(M2- 1)3/2- 

±A/2(1 + Cp(l)y 





±M 




±(A'/2 _ 1) 




_ - 1 _ 



M3(1 + cp(l)) - M(M2 - 1) + A/(A/2 - 1) il/3(l + cp(l)) 



(A/2 - 1)3/2 



(A/2 - 1)3/2 



SO by (|72|) the it-ficlds are genuinely nonlinear at ?7. If we choose e > 
sufficiently small, then they are genuinely nonlinear for all values of U with 
\\U-U\\<e. 

Consider the upper left quadrant, y > Q > x, Here A_ is relevant. In 
increasing x direction with fixed y, corresponding to decreasing ^, the change of 
i7 in a simple wave is given by — r_: density increases, velocity turns down and 
decreases (same effect as the A_ -shocks). This is a compression wave. It can be 
approximated as the limit of an increasingly fine fan of weakening shocks. 

In the upper right quadrant x,y > 0, \+ is important. In increasing x 
direction with fixed y, corresponding to decreasing ^, the change of J7 in a 
simple wave is given by —?'_(.: density decreases, velocity turns downwards and 
mcreases (opposite to the behaviour of A+-shocks). This is an expansion wave 
(also known as Prandtl-Meyer wave). 



18.2 Summary 

All results combined, we have the following description of steady and self-similar 
Euler flows U that are sufficiently L°°-close to a constant background state 
U = {p,Mc,0) with Mach number A/ > 1 (supersonic), defining Mach angle 
p, = arcsin (see Figure [T|) : 

1. they are necessarily BV, 

2. they are constant outside six narrow sectors whose center lines are (1 : 0), 
{cosp : sin/i), (cos// : — sin/t). 
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3. in the (1 : 0) forward and backwards sectors U is constant on each side 
of a single contact discontinuity (which may vanish) , 

4. in the forward (cos/z : ±sin^) sectors U is constant on each side of a 
single shock or single rarefaction wave (which may vanish) , 

5. in the backward (cos/x : ±sin^) sectors U can have an infinite or any 
finite number of shocks and compression waves, but 

5a. two consecutive compression waves with a gap are not possible, and 
5b. the shock set (on the unit circle) is discrete, with each shock having con- 
stant neighbourhoods on each side whose size is lower-bounded proportionally 
to the shock strength. 

It does not seem possible to improve these results without making additional 
assumptions. Examples with infinitely many consecutive shocks, or shocks in- 
terspersed with compression waves, or compression waves ending in a point that 
is a limit point of shocks, can be constructed. 

19 Appendix 

19.1 Saltus functions 

Definition 13. A saltus function f : D ^ IR"' (D C IR) has the form 

/(^) = ^ an+ ^ bn 

Xn<X X„<X 

where (x„) is a sequence and X) absolutely converging series. 

19.2 Versions 

Lemma 14. Let ft C IR" measurable nonempty, K C IR™ compact, U E L°°{il) 
so that U{x) G K for a.e. x G fl, g : K — > IR*^ and g : £7 — > IR'^ continuous. If 

g{U{x)) < g{x) for a.e. x eD,, (74) 

(meaning gi{U {x)) < gi{x) for all i, where g = {gi, ...,gk), 9 = {gi, -.gk)), then 
we can find a version U of U , with values in K everywhere, so that 

90) < g for all x e n. (75) 

Proof. We immediately modify U, on a set of measure 0, to have values in K 
everywhere. 

Let E = {x\ g{U{x)) < g{x)}. Then CE has measure zero, so every x € CE 
is the limit of a sequence (a;„) in E. {U{xn)) C K which is compact, so we may 
choose a subsequence {x'„) so that ([/(x^)) converges as well. Define U{x) := 
lun{U{x'„)) G K (we use one subsequence for each x, as the limit for others may 
be different of course) . Now 

g{U{x)) ^ g{U{xn)) < g{xn) -> g{x), 
so ((75)) is satisfied. □ 
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